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INTRODUCTION 


This report documents the background information in establishing 
several proposed atmospheric turbulence models for use on motion based 
aircraft simulators. Simulated time histories of aircraft motion in a 
turbulence environment are required in a variety of engineering applica- 
tions, and their use appears to be Increasing as more intricate and 
sophisticated design studies are attempted. As an example, the use of 
flight simulators for the study of airplane handling qualities and ride 
quality has proven to be more valuable when disturbances in the form of 
artificially simulated turbulence are introduced into the system. Several 
methods have been used to generate turbulence signals; each one aimed at 
realising the actual atmosphere as closely as possible. A realistic 
representation of turbulence becomes especially important in the simula- 
tion of future aircraft with high sensitivity to turbulence, as even 
light to moderate turbulence may seriously degrade their controllability 
and ride quality. The iow altitude atmospheric turbulence critically 
effects the evaluation of vehicle handling qualities, pilot work load, 
ride quality, and other design factors. Several emperical studies 
(Refs. I, 2, 3) have shown that low altitude clear air atmospheric tur- 
bulence is only locally isotropic i.e. isotropic over a finite range of 
wavelengths. In order to account for the anisotropy of typical low 
altiTude clean air turbulence, the rms velocities of the gust field are 
randomly varied in the proposed gust model. The proposed model, in 
addition to varying turbulence intensity (rms velocity), varies the 
atmospheric turbulence scale length. The scale lengths predicted by 
either the Von Karman or the Drv'uen models (Ref. 5) are large compared 
to real atmospheric turbulence. The scale length distribution Is, there- 
fore, modified to achieve compatibility with real atmospheric turbulence. 

With a suitable combination of scale length and intensity distribu- 
tion, the proposed model wilt simulate various atmospheric conditions 
characterized by altitude, stability, and ter»-ain. This new model is 
mechanized to be included in a flight simulator experiment in order ro 
determine to what extent the pilots are sensitive to changes in atmos- 
pheric conditions and the realism of the model. 



The fol lowing sections of this report describe the proposed turbu- 
lence model and the flight simulator experiment in detail. Briefly the 
sections consist of: 

(a) Literature Survey: Since atmospheric turbulence is 3 stochastic 

process, a review of probability and statistics is included in this sec- 
tion. In addition the statistical properties of real atmospheric turbu- 
lence is discussed. 

(b) Presently Used Techniques: Presently used techniques are inum- 

erated and the need for a new model is demonstrated. 


(c) Proposed Model: The proposed turbulence mode! is discussed 
and the theoretical results compared with real atmospheric turbulence 
demonstrating an improved representation of atmospheric turbulence. 


(d) Simulation Details: 
flight simulator experiment in 
of the various turbulence mode 
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LITERATURE SURVEY 


In this section a review of probability and statistics is followed 
by a summary of the statistical properties of atmospheric turbulence. 

Review of Probability and Statistics (Refs, 8^ 15) ; 

Stat lonarl ty : A random process is stationary if its statistical 

properties are not dependent on the time of their measurement. One could, 
for example, collect an infinite number of time histories, called an 
ensemble, which are representative of the process. If one takes an avei — 
age across the ensemble, and if these averages are not a function of time, 
the process is stationary. 

Homogeneity: A random process is homogeneous if its statistical 

properties are independent of position. 

Ergodicity: In turbulence measurements It is impossible to obtain an 

ensemble from atmospheric measurements. Thus it is necessary to use time 
averages to get statistical information, if such a time average yields 
the same statistical properties as the ensemble average the process is 
caJ led ergodic. 

Mean Value: The mean value of a random variable, u, of an ergodic 

random process is given by 

u = Lim /^-j- u(t )dt 

T-><» “ 

In practice the limit is not required and u can be approximated by 

u ~ Y !q u(t)dt, for T large. 

This approximate representation is especially useful for processes such 
as turbulence. However, the time interval T must be large enough so 
that the average approaches the asymptotic value one would obtain for 
a stationary process. 

Variance: The variance of u is defined as 

= Urn ^ C(u(t) - u)^ldt. 
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as before in practical applications the variance can be approximated by 


“S = T Cu(t) - 


for sufficiently large T. 

Standard Deviation: The standard deviation is defined as the 

sqiare root of the variance. 

th 

Normalized Centra! Moment: The n normalized central moment, M of 

a random process, u(t), is 


M = Lim 
n 2T 


/ItC' 


uCt) - U -|2 


■Tdt, n = 1,2,3, 


which can be approximated by 




— ]*dt, n = 1,2,3, 


Cumulative Probability Distribution: The cumulative probability 

distribution of u(t), P^(x) is defined as the probability that u £ x. 

Probability Density Distribution: The probability density distribu- 

tion of u(t), P^j(x) is defined as the probability that; x < u < x + dx. 

Gaussian Probability Density Distribution: If a random variable, 

u(t), is Gaussian distributed its probability density is given by 


P (x) = -J — 
u 


expC“ 


I (x - u) 
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Rayfeigh Distribution: Another probability density of interest is 

the Rayleigh Distribution defined as follows: 


P(x) = exp(“ “ x^/c^) 
c^ 2 

where c^ is one half the expected value of the random variable x or 


c 


2 
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Cross Correlation Function: The cross correlation function of two 

random processes u(t;, w(t) is defined as 




I 

2T 


R (t) = Urn u(t)w(t + T)dt 

UW J -T 


correlations are the measures of the predictability of a signal at some 
future time (t + t) based on the knowledge of a signal at time t. 

Autocorrelation Function: The autocorrelation function is a special 

case of the cross correlation function defined above in which w(t) = u(t), 
such that, 


R (t) = 
uu 


L'm 2T -^-T 




u(t)v(t + T)dt, 


Integral Scale Length: A statistical parameter of special importance 

in atmospheric turbulence is the integral scale length, 


L 

u 



r 

J .0 


R (x}dT. 
uu 


where Uo is the reference steady state flight speed. 

Cross Spectral Density: The cross spectral density of two random 

processes u(t) and w(t) Is defined as the Fourier transform of their cross 
correlation. 


$ Cf) = / R (T)expC-i2TTfT)dT, 
where f Is frequency. 

Power Spectral Density: The power spectral density, PSD, of a 

random process is the Fourier transform of its autocorrelation function, 
or 


(1) (f) = R (T)Exp( i2'irfT)dT. 

The PSD can be Interpreted physically bs the average contribution to the 
variable from the frequency component f. Thus, 
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d>2 = /“ A (f)df. 

White Noiset White noise Is a random process for which the PSD is 
a constant Independent of frequency. That Is, (f)o(f) = k. 

Properties of Atmospheric Turbulence : AtriKDspheric turbulence Simula-' 

tion studies, in general, begin with the study of the real atmosphere. In 
Refs. I, 2, 3, 4, atmospheric data have been reported characterizing 
various atmospheric conditions In the form of terrain, stability, altitude, 
temperature, time, season and geographic location. This data has been 
suitably modified to establish a basis of comparison for the simulated 
turbulence field. 

The fol lowing criteria are used for the basis of comparison: 

Output Statistics: 

Mean 

Standard Deviation 
Probability Distribution: 

Cumulative Probability 
Probab i 1 1 ty Dens i ty 
Fourth and Sixth Moment 
Patchiness of the Field; 

Power Spectral Density: 

Element of Surprise. 

Each of these properties will be discussed from the standpoint of real 
atmospheric turbulence. 

Mean: u (mean velocity of ' ;^Itudinal turbulent gust component) 
u /q u(t)dt where T I::, large. 

The mean of a real atmospheric turbulence field is 0 ± 0,1 ft/sec. 
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standard Deviation: (standard deviation of longitudinal turbulence 

gust component) of a random process is approximated by: 

- u)^dt 

u T ■'0 

The standard deviation of the velocity field for low altitude clear air 
turbulence Is 3.0 ± 1.51 ft/sec. 

Typical values are also listed in Ref. 5 

For clear air turbulence: 

0=2 ft/sec for light turbulence, 
u ^ 

0=4 ft/sec for moderate turbulence, 
u 

o^ = 6 ft/sec for heavy turbulence, and 

o, =0=0 =21 ft/sec for thunderstorms, 

u V w 

Probability Distribution: (Reference 8 ). 

A typical cumulative probability distribution of atmospheric turbu- 
lence indicates a departure from a Gaussian process, showing increased 
probability of both large and small gusts. The same is true of the prob- 
ability density function. 

A table beiow compares the normalized fourth and sixth moments of 
atmospheric turbulence to a Gaussian process. 



M 4 

Ms 

Atmosphere 

3.5 

21 .7 

Gaussian Process 

3.0 

15.0 
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Patchiness: it is known that turbulence has a patchy structure, which 

seems to occur in bursts of relatively intense motion separated by areas 
of relative ca Im. 

Power Spectral Density (PSD): Analysis of data for about 40 turbulence 

fields (Ref. I) characterizing various atmospheric conditions shows that at 
high frequencies the spectral density varies as A constant horizon- 

tal asymptote seems to fit best at lower frequencies for the longitudinal 
and lateral components. 

Element of Surprise: More often than not, real atmospheric turbulence, 

when enccuntered, presents an element of surprise. It is not easy to 
formulate a model of this phenomenon in terms applicable to flight simu- 
lator work. It seems that a measurement of "sudden Jump" In the velocity 
field cari be used as a possible criterion to describe this phenomenon. 

The relative frequency of "sudden jump" of atmospheric turbulence can be 
compared to the simulated model. Changes in aircraft orientation angles 
can also be used to measure this phenomenon. 

Presently Used Simulation Techniques: In the preceding pages of this 

report we have discussed the slatlstical properties of atmospheric 
turbulence which are to be modeled by a realistic simulation. In this 
section several presently used simulation iachniques are discussed from 
the standpoint of their statistical realism and suitability for use in flight 
simulators. 

Measured Turbulence Field: Flight record ‘ ngs of atmospheric 

turbulence is perhaps the most obvious method of producing a realistic 
simulation. There can be little argument as to whether or not these time 
histories are an accurate and re-jlistic representation. However, it is 
difficult to adjust the measured time histories to allow for conditions 
other than those for which it was recorded. No allowances can be made 
for changes of altitude or different atmospheric conditions. Another 
serious drawback is that the recorded time histories are fixed in length. 
Extended run times, therefore, cannot be accommodated without repitition. 

From the simulation point of view the pilots tend to recognize some of the 
characteristics of the turbulence field and develop an intuition for pre- 
dicting the field. This defeats the purpose of an artificially simulated 
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turbulence field, which is to provide unpredictable external disturbances. 
It can, therefore, be concluded that flight recordings of atmospheric tur- 
bulence are not suitable for the simulation of typical turbulence. 

Sum of Sine Waves; Reference 8 describes this method in summary 
form. This technique Involves superimposing several sinusoidal waves 
of different frequencies and amplitudes. The resultant Is used to repre- 
sent time histories of turbulence. One obvious disadvantage of this 
method Is that it contains only a finite range of frequencies whereas 
actual atmospheric turbulence consists of an infinite number of frequency 
components. 

Results of this simulation are not available but the model can justi- 
fiably be discarded on the basis of its inadequacy in matching the fre- 
quency content. 

Method of Orthogonal Functions; In this method (Ref. 4) the 
recorded time histories of turbulence are decomposed into eigenfunctions 
of a covariance matrix. The probabilistic structure of the eigenfunction, 
and the coefficients of each of the time histories are studied. Simulated 
time histories are then regenerated by suitably modifying the distribution 
of the coefficients. The available preliminary results show that this 
technique adequately models the frequency contents and also presents an 
element of surprise. However, this model falls to show a patchy non- 
Gaussian characteristic which Is typical of the real atmosphere. In 
addition to the mathematical complexity of the technique, its application 
is limited since recorded time histories are needed, 

Gaussian Turbulence Model; The classical method, most widely used 
1 3r turbulence simulation, is the linearly filtered white noise technique. 
Here the turbulence gust field is produced by passing white noise through 
a linear filter as shown in Figure la. The resultant signal is shaped so 
that the power spectrum and rms intensities match those of real turbulence. 
A Dryden or Von Karman form (Ref. 5) are normally used to model the power 
spectrum. This model is remarkably easy to implement and can be adjusted 
for any general power spectrum. However, this model too falls short of 
reproducing the non-Gauss Ian patchy nature of real turbulence. Figure ib 
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Gust 

■> Velocity 
Component 




40 Sec . Samp I e 
Figure ib 


Derivative of Vertical Velocity 




compares the artificially simulated gust field using this Gaussian model 
(with a Dryden spectrum) and real atmospheric turbulence. It may be 
observed that the intensity for the Gaussian model is nearly constant 
whereas measured ("rea!") turbulence exhibits a patchy nature or intensity 
bursts. Test pilots, when exposed to this model in a flight simulator, 
rated the realism fair to poor (Ref. II). 

Non-Gaussian Turbulence Model: References 7, 8, It present a non- 

Gaussian turbulence model. Time histories are generated by muitiplying 
two independence random variables, one to represent the turbulence within 
a patch and the other to represent the variation of intensity with time. 
Figure 2 shows two independent (Gaussian white noise generators and linear 
filters, which produce Gaussian random variables, a(t) and b(t). These 
variables are then multiplied to produce gust time histories. 

The non-Gaussian model proposed in Ref. 8, a modification of the 
above, Is shown in Figure 3. Here a(t), b(t), and d(t) are Independent 
Gaussian processes. The process c(t) fs generated by multiplying a(t) 
and b(t). The resultant process, c(t), a modified Bessel process, ‘3 
summed with d(t) to form the output, u(t). The most remarkable achieve- 
ment of this model is that the patchy characteristic and several statis- 
tical parameters of the simulated turbulence field can be varied simultan- 
eously by varying the standard deviation ratio (R = hfowever, 

when R Is varied to achieve one set of statistical properties, several 
other statistical parameters of interest do not match real turbulence. 

In addition, due to the mathematical complexity, the mechanization of 
this mode! on a flight simulator is complicated and expensive (Ref. 8). 

It can be observed from the review of presently used simulation tech- 
niques that there is a need for a new model which adequately matches real 
atmospheric turbulence and is simple to implement in flight simulator 
studies. None of the preceeding models have the flexibility of simulating 
various atmospheric conditions characterized by altitude, stability, and 
terrain. It Is, therefore, necessary to introduce a new turbu'ence model 
which is realistic and can flexibly accommodate changes in atinospheric 
conditions and be easily implemented in flight simulator studies. 
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Figure 2 

Non-Gaussian Turbulence Model (Reference \) 
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PROPOSED GUST MODELS 


Of tha simulation techniques described, the Gaussian turbulence 
model is the simplest to Implement and least expensive computat Iona 1 1 v . 

The proposed models, modifications of the Gaussian simulation technique, 
retain the simplicity of the Gaussian technique while modeling the 
characteristics of real atmospheric turbulence. In this report three 
basic models are proposed. 

1) Modified Gaussian Model 

2) Rayleigh Model 

3) UVA Turbulence Mode! 

Modified Gaussian Model: A block diagram of the Modified Gaussian 

Mode! Is presented In Figure 4a. Gaussian white noise, (f)o. Is passed 
through a linear filter, G(S), whose power spectrum is given by a Dryden 
model, in order to avoid computational complexity the Dryden form is 
selected In this report over the Von Karman form. The linear filter, G(S), 
is modified to include random variations of rms intensity. The random 
number generated by A is passed through a distribution modifier to gen- 
erate rms intensities. Time histories are then generated by passing 
Gaussian white noise, (J>o, through the linear filter modified by the 
distribution modifier. 

The patchy nature of atmospheric turbulence suggests that the field 
is composed of two componenhS: One to represent variation of intensity 

within a patch and the ether to represent variation of intensity with 
time. The distribution modifier in +his model, essentially, represents 
the variation of intensity with time. The level of turbulence within 
each patch is controlled by the magnitude of the rms intensity. 

The Distribution Modifier Is the probability density function of 
the rms intensity. Analysis of several sets of atmospheric data characterized 
by various atmospheric conditions show that a truncated Gaussian distribu- 
tion best fits the probability density of rms intensity (Reference 1). In 
this report two sets of data characterized fay terrain, altitude and 
stability are derived and presented in Table la along with the atmospheric 
conditions. Throughout this report the gust field generated by these two 
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Figure 4a 

Modified Gaussian Turbulence Simuiation 
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Table I A 


DISTRIBUTION MODIFIERS 


RMS DISTRIBUTION 
MODIFIER 
CASE It 

0 ft/ sec 
u 

cr ft/sec 

V 

a ft/ sec 
w 

MEAN 

3. 1 

3.2 
2.8 

VARIANCE 

1.2 
1 .2 
0.9 

RMS DISTRIBUTION 

a ft/sec 
u 

3.2 

0.8 

MODIFIER 

a ft/sec 

V 

3.5 

1.0 

CASE f ! 1 

0 ft/sec 
w 

4. 1 

0.0 1 


CASE 1 I : 


CASE 111: 


ALTITUDE; 

STABILITY; 

TERRAIN: 

ALTITUDE; 

STABILITY: 

TERRAIN: 


250 Ft. 

UNSTABLE 

PLAINS 

750 Ft. 

UNSTABLE 

MOUNTAINS 


17 






Distribution Modifiers will be referred to as Case fl and JM. (Case I 
is a Gaussian model.) 

Rayleigh Model: The Rayleigh model Is derived from the previous 

model by replacing the Distribution Modifier by a Rayleigh probability 
density function. The Rayleigh probability density function for rms 
vertical turbulence intensity, a. . is given by 

W 

p(a^) = -^Exp(- ia^c^) 
c 

where c^ Is one half the expected value of a^. 

w 

Rrom the Dryden spectrum models of real atmospheric turbulence the 
value of c has been estimated in Ref. 5 to . 2.3 ft/sec. 

The rms intensities of the longitudinal, u, and the lateral, v, gust 
components are obtained from the relation: 



ILL 
U V w 

The scale lengths are given by: 

L = L = L h > 1750 ft. 

u V w 

ya 

= 145 h h < 1750 ft. 

L = h. 
w 

This model will be referred to as Case IV. 

UVA Turbulence Model: The UVA Turbulence Model includes in addition 

> 

to the rms Distribution Modifier, a scale length modifier. A block 
diagram of this model is presented In Figure 4b. In addition to con- 
trolling the patchiness of the turbulence field, the time variations of 
scale length achieves numerical compatibility with the real atmosphere 
and further randomizes the simulation. 
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Figure 4b 

UVA Gust Model Turbulence Simulation 
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The Scale Length Distribution Modifier is derived from data collected 
In the LO-LO-CAT program (Ref. I) for various combinations of altitude, 
terrain, stability, temperature and geographic location. Figure 5 shows 
the Gaussian distribution of scale length derived from the Ref. I. Of 
the several atmospheric conditions analyzed. Table lb presents the scale 
length modifier along with the corresponding rms modifier. This model 
will be tested for two separate atmospheric conditions characterized by 
altitude, terrain and stability. These two cases will be referred to as 
Cases V and VI. 

Attempts are also being made to fit an accurate non-Gaussfan distri- 
bution for both scale length and rms distribution modifiers. 
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Scale Length - Ft 
Figure 5 

Scale Length Distribution 
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Table lb 


DISTRIBUTION MODIFIERS 



1 

Case 5 

Case 6 



MEAN 

VARIANCE 1 

MEAN 

VARIANCE 

RMS Distribution 

cj ft/sec 

3.! 

1 .2 i 

3.2 

0.8 

Modifier 

u 

a ft/sec 

V 

3.2 

[ .2 

3.5 

1 .0 


a ft /sec 
w 

2.8 

0.9 

4. i 

0.9 

Scale Length 

L ft 

415.0 

1 10. 0 


1 16.56 

Di stri but ion 
Modi f ier 

u 

L ft 

V 

325.0 

86,65 

wm 

126.64 


L ft 
w 

535.0 

83, 1 1 

I 425.0 

i 

132.98 


Case 5 

Altitude: 250 ft 

Stability: Unstable 

Terrain: Plains 

Case 6 

Altitude: 750 ft 

Stability: Unstable 

Terrain: Mountains 
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RESULTS 


In this section results obtained by statistical analysis of the gust 
velocity components for each of the five models will be discussed In view 
of the properties of real atmospheric turbulence. The statistical results 
have been obtained in the form of: 

1) Mean and standard deviations. 

2) Normalized fourth and sixth moments. 

3) Probability density functions, 

4) Power spectral densities. 

5) Patchiness. 

6) Relative frequency of element of surprise. 

Table 2 tabulates the mean and standard deviation of gust components 
of each of the six cases to be simulated. It can be observed that the 
standard deviation varies from 2.6 to 5.2 ft/sec which is typical of low 
altitude clear air Clight to moderate) turbulence. 

Fourth and sixth moment characteristics are tabulated in Table 3. 

Within the limits of experimental error these characteristics are in 
fairly good agreement with the real atmospheric data obtained in the 
LO-LO-CAT experiments {Ref. I). 

Since the cumulative probability and the probability density function 
essentially contains identical information, only probability density 
functions will be analyzed. Figures 6 to 10 are plots of probability 
density functions for the simulated cases, in order to compare with 
atmospheric turoulence, a Gaussian distribution is plotted on the same 
scale. It has been established (Ref. 8) that real atmospheric turbulence 
exhibits a higher probability of both smaller and larger gust velocities 
than the Gaussian distribution. A careful study of the probability density of 
thesimulated field reveals a higher probability of larger gust velocities 
compared to a Gaussian distribution, however, the distributions do not show 
higher probability of lower gust velocities. 
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Table 2 


MEAN AND STANDARD DEVIATION OF 
GUST COMPONENTS FOR PROPOSED SIMULATION 
(over a 2 mile segment) 


CASE 

OUTPUT STATISTICS 

GUST COMPONENT 

MODEL 



UG 

VG 

WG 


1 

MEAN 

0.086 

0.063 

“0.031 

GAUSSIAN 


STANDARD DEVIATION 

3,97 

3.90 

4.43 


2 

MEAN 

0.83 

“0.32 

-0.15 

MODIFIED 


STANDARD DEVIATION 

3.9 

3.5 

2.6 

GAUSSIAN 

3 

MEAN 

0.88 

“0.40 

0.06 

MOD 1 F 1 ED 


STANDARD DEVIATION 

3.9 

3.9 

3.8 

1 

GAUSSIAN 

4 

MEAN 

-0.36564 

-0. 1663 

-0.2256 

RAYLEIGH 


STANDARD DEVIATION 

5.1998 

4.8406 

4.4880 


5 

MEAN 

0.27621 

-0.3604 

-0.2080 

U, Va. 


STANDARD DEVIATION 

3.6644 

3.5022 

2,6756 

MODEL 

6 

MEAN 

0.2010 

“0. 1017 

“0.3261 

U. Va, 


STANDARD DEVIATION 

3,5680 

3.9064 

3,8100 

MODEL 
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Table 3 


FOURTH AND SIXTH MOMENT DATA OF 
REAL AND SIMULATED TURBULENCE FIELDS 
Cover a 2 mile segment) 


RFAL ATM. 

FOURTH 

SIXTH 

3.5 
21 .7 

3.5 

21.7 

3.5 
21 .7 




GUST ) 

/ELOCITY C0^ 

^PONENTS 



MOMENT 

UG 

VG 

WG 

MODEL 


FOURTH 

3.0 

3.86 

3.00 


1 

3. 




GAUSSIAN 


SIXTH 

15.0 

15.60 

15.00 


2 

FOURTH 

5.864 

5.546 

3.176 

MODIFIED 


SIXTH 

61.076 

22.323 

16.846 

GAUSSIAN 

3 

FOURTH 

5.129 

3.224 

2.853 

MODIFIED 


SIXTH 

46.698 

18.976 ^ 

1 

1 J .985 

GAUSSIAN 

4 

FOURTH 

3.738 

3.200 

3.3085 

RAYLEIGH 


SIXTH 

21.726 ; 

18.075 

19.887 ' 

1 


5 

FOURTH 

3.494 

3.236 

3.467 

1 

U. Va. 


S 1 XTH 

20,838 

16.012 

21.803 ; 



FOURTH 

3.065 

3.145 

3.856 ' 


^ i 





U. Va. 


SIXTH 

14.006 

16.12 

21.475 
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Probabi I Ity 








0 


3.75 


7.75 


I I .75 


Gust Velocity Component 
w - Ft/Sec 

g 

Figure 9 

Probability Density (Case V) 
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Gaussian Process 
S i mu I ated Field 



Gust Velocity Component 
w - Ft/ Sec 

g 

Figure 10 

Probability Density (Case VI) 
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Power spectral densities of the simulated turbulence fields are 
presented In Figures II to 15. The high frequency contents are compared 
with a line of slope -2. It is observed that the power spectrum at high 
frequency varies as (-2 logarithmic slope). At low frequencies a 
constant asymptote fits best. The power spectrum in the entire frequency 
range within the limit of expe'-imenta I error, is in fairly good agreement 
with the assumed Dryden form (Figure 16). 

The patchiness of each of the cases is presented in Figures 17 to 19. 
The derivative of the vertical gust component is plotted illustrating a 
varying intensity of patchiness. Case IV presents patchy characteristics 
which closely matches real atmospheric turbulence. 

Element of surprise is tabulated in Table 4 . The presently used 
criterion, "sucden jump in velocity field,” does not adequately quantize 
this phenomenon. For future work an alternative criterion namely, changes 
in aircraft orientation angles, may be used to measure this phenomenon 
quantitatively. 
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Case VI 



Figure 16 

Patchiness of Simulated Field 
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Figure 19 


Block Diagram of Motion Based Simulator (Reference 15) 










Table 4 


FREQUENCY OF ELEMENT OF SURPRISE 
OF SIMULATED FIELD 


Case 

FREQUENCY OF ELEMENT OF SURPRISE % 


UG 

VG 

WG 

2 

0.0666 

0.0666 

0.0 

3 

0.0333 

0.0 

0.0 

4 

0.0 

0.0 

0.233 

5 

0.0533 

0.0333 

0.0 

6 

1 

0.2666 

0.400 

0.0 
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SIMULATION DFTAiLS 


This section describes the flight simulator experiment, including Jetaiis 
of the aircraft simulated, the flight simulator and the pilot task perform- 
ance. The purpose of the experiment is to determine to what extent pilots 
are sensitive to changes In atmospheric conditions and realism of the 
simulation. 

The aircraft simulated in this study is the Canadian de HavMland 
DHC-6 Twin Otter. This particular aircraft is chosen because it is a 
typical STOL aircraft and Its flying characteristics are well known. 

In addition, there are many pilots available with flying experience in 
the Twin Otter to valldidate the simulation. Table 5 tabulates typical 
aircraft parameters (Ref. 15). 

The simulator used in this study is the six degree of freedom visual 
motion simulator (rms) at the NASA Langley Research Center. This is a 
motion based simCitator with the basic interior and instrumentation of a 
jet transport cockpit. Figure 19 presents a block diagram of the simu- 
lator reproduced from Ref. 15. 
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Table 5 


AIRCRAFT PARAMETER (Ref. 15) 


w = t 1500 lb 

u = 256.67 ft/sec 
o 

= 0.045 

h* = 0.2 

1^ = J6900 siug-ft^ 
[ = 27600 slug-ft^ 


1^ = 40600 slug-ft^ 
F = 1400 slug-ft^ 

c = 6.5 ft 

b = 65 ft 
s = 420 ft^ 


Pf<-0T POSITION WITH RESPECT TO C.G. 
X = 8.8 ft y = “1 .6 ft 


2 
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DESCRIPTION OF THE EXPERIMENT 


The tests will be conducted for six cases composed of Gaussian 
turbulence and five models. Table 6 presents preliminary planning of 
the test runs along with the duration and purpose of each run. The 
pilots will fly each of the six runs In random order for ten minutes 
achieving a composite flight task. Initially the pilot. In cruise, will be 
required to achieve a constant altitude tracking task. After each run 
the pilot will be given a flight questionnaire (see Appendix A) to deter- 
mine sensitivity to the various models and realism of the simulation. The 
entire experiment will be repeated with a higher level of turbulence 
intensity. 



Table 6 




PRELIMINARY PLANNING 

OF SIMULATOR 

EXPERIMENT 


SUB 

RMS DISTRIBUTION 

SCALE LENGTH 

ATMOSPHERIC 

PURPOSE OF RUN 

DURATION OF RUN 

MODEL 

MODIFIER 

DISTRIBUTION MODIFIER 

CONDITION 



1 

DETERMINISTIC 

DETERMINISTIC 

AVERAGE 

TO EXPOSE THE 

10 Min 


ff = 4.0 
u 

L|j = 1317.41 


PILOT TO GAUSSIAN 
MODEL 



0 = 4.0 

V 

L^ = 1317.41 





a = 4.5 

VJ 

L = 750 
w 




2 

GAUSSIAN 

DETERMINISTIC 

250 FT, 

TO EXPOSE THE 

10 Min 


DISTRIBUTION 


UNSTABLE 

P 1 LOT TO EXTREME 



MEAN ST. DEVIATION 

L = 913.44 
u 

PLAINS 

ATMOSPHERIC 

CONDITIONS 



U: 3.1 1.2 

L = 913.44 

V 





V: 3.2 1.2 

L = 250 
w 





W: 2.8 0.9 




3 

GAUSS 1 AN 


750 FT, 

TO EXPOSE THE 



MEAN ST. DEVIATION 
U: 3.2 0.8 

L = 1517.41 
u 

UNSTABLE 

MOUNTAINS 

PILOT TO EXTREME 

ATMOSPHERIC 

CONDITIONS 

1 0 Mi n 


V: 3.5 1.0 

L = 1317.41 

V 





W; 4.1 0.9 

L = 750 





w 



4 


RAYLEIGH 


DETERMINISTIC 


C = 2.3 L = 913.44 

u 

= 913.44 

L = 750 
w 

5 GAUSSIAN GAUSSIAN 

MEAN ST. DEVIATION MEAN ST. DEVIATION 
U 3. I i.2 U 415 1 10.0 

V 3.2 1,2 V 325 86.65 

W 2.8 0.9 W 335 83.11 

6 GAUSSIAN GAUSSIAN 


UNKNOWN 


TO PROVIDE FEEL 
FOR PATH I NESS 


10 Mi n 


250 FT, TO EXPOSE THE 

UNSTABLE P I LOT TO EXTREME 

PLAINS ATMOSPHERIC 

CONDITIONS 10 Min 


250 FT, TO EXPOSE THE 

UNSTABLE PILOT TO EXTREME 

MOUNTAIN ATMOSPHERIC 

CONDITIONS 10 Min 
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APPENDIX A 


FLIGHT QUESTIONNAIRE: 

FI ight Number Date 

Pilot: 

1. Turbulence Intensity: 

Light Moderate Severe Extreme 

2. Realism of Turbulence: 

Very Good Good Fair Poor Very Poor 

3. Correctness of Relative Amplitude of Disturbances: 

Not Enough About Right Too Much No Comments 

Roll 

Pitch 

Yaw 

Heave ___ 

Side Force 

4. Patchy Characteristics (Variation of Intensity Bursts) 

Much Too Continuous A Little Too Continuous About Right 

A Little Too Patchy No Comments 

5. Frequency Contents of Turbulence: 

Not Enough About Right Too Much No Comments 

Low FRQ: 

High FRQ: 

6. Element of Surprise in the Simulated Turbulence Field: 

a. Quite Often Sometimes Never 

b . Rea I i sm of 6a : 

Very Good^_ Good Fair Poor Very Poor 

7. Atmospheric Conditions: 

a. Altitude: 0 - 1,000 Ft 1,000 - 10,000 Ft 

Over 10,000 Ft ^ Unable to Judge 
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8. Pilot Estimate of the Work Load: 

Very Easy Easy Average Difficult Very Difficult 

9. Pilot Estimate of Task Performance: {Integral Squared Error for ILS 

Tracking Task) 

Very Good Good ___ Average Poor Very Poor 

10- Realism of This Model Compared to Previously Flown Model: 

Very Good Good About the Same Poor Very Poor 

11- Did You Observe a Repetitive Pattern In the Turbulence Field? 

Yes No 

12. Cooper-Harper Rating; 

13. Additional Comments About Realism of Turbulence and Aircraft Simulation: 
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PILOT EXPERIENCE: 


1 . Name Date 

2. What Type of Flying Experience Have You Had? 

Military Civil 

3. Main Types of Aircrafts Flown: 


4. Total Number of Hours Flown: 

5 , Hours of Instrument Flying: 

G. Hours in Simulators: 

7. Hours in VMS: | 

I 

J 

8. Hours in Twin Otter; 

9. a. Estimate the % of Time Flown in Turbulence: 

b. Of This Time What % Was Flown In 

j 

Light Turbulence Moderate Turbulence Severe Turbulence Extreme Turbulence ] 


10. What Characteristic of Turbulence I nterf eres Most with Your Ability to 
Control the Aircraft? 


11. Describe the Most Critical Case of Turbulence Encountered During Your 
Flying Experience: 

a. Day Night 

b. Terrain: Altitude: 

c. Atmospheric Stability: 

Stable Neutral Unstable Unable to Judge 

d. What Was fhe Task You Were Attempting Before Turbulence Was Encountered: 

(e.g. ILS Approach, Cruise, etc.) 

e. Any Additionaf Comments: 

•i 
f 

r 
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